Effective degrees of freedom during the radiation era 
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We update the curves of the effective degrees of freedom for energy density g*(T) and for the 
entropy density g*s(T) during the era of radiation domination in the Universe. We find that a plain 
count of effective degrees of freedom sets an upper limit to the temperature of the quark-hadron 
transition at T c < 235 MeV for the energy density and T c < 245 MeV for the entropy density. 
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As is well known, the thermal history of the Universe 
during the radiation era depends on the effective degrees 
of freedom g e f / (T) or gr* (T) contributed by the particles 
in thermal equilibrium at temperature T. The function 
<7*(T) was first computed over the eV to TeV range in 
1981 Q and later updated in 1988 and 1990 0, and 
has often been quoted without updating since then (e.g. 
4], 5j). Over time our knowledge of the properties of 
the particles determining this function have undergone 
various modifications, particularly at high temperatures 
due to improved values for the top quark mass, improved 
limits for the Higgs boson, and improved detail on the 
hadron spectrum ||j . This calls for an updating of g* (T) . 

Recall that the light nuclei arc synthesized at temper- 
atures T <C m e when all particles have decoupled and 
g» receives contributions only from the three relativistic 
neutrino species and the photon. Each neutrino species 
then contributes 7/4 degrees of freedom, corrected for the 
evolution of the neutrino temperature through the period 
of e annihilation, and the photon contributes g 7 = 2. 
The function g* (T) then has the low-energy limit 



The total energy density of all species in equilibrium 
at photon temperature T can be written Q 



9 = m g * T ■ 



(2) 



Summing over all species i with mass m^, the effective 
degrees of freedom g* for the energy density is given Q 
by 



9* = 
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(3) 



Here Ti is the temperature of species i if it has a thermal 
distribution with a different temperature than photons 
(e.g. neutrinos), Xi = rrn/Ti, and the ± sign implies + 
for fermions and - for bosons. A difference from Kolb 
and Turner Q is that we neglect the chemical potential 
which one can do for a thermal radiation background. 
The entropy density is defined by 



g* = g~/ + 3 



4/3 



3.363 . 



(1) 



At temperatures above this low-energy limit but in the 
'confined' phase below the quark-hadron phase transition 
T c , the hadrons contribute by the pions, g w — 3, and by 
the low-energy tails of the heavier mesons and baryons. 
At all temperatures the three charged leptons contribute 
by gi = 3 x |, and the vector bosons by gw± = 6 
and gz« = 3. In the 'deconfined' phase above T c the 
hadron contributions are replaced by the six quarks, each 
of which contribute g q = 12, and by the gluons, g g = 16, 
and the Higgs boson, g H o — 1. We take their masses 
from the Review of Particle Physics which, for the 
u, d and s quarks implies current-quark masses. 
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(4) 



where g^s is the effective degrees of freedom . The con- 
tribution from all species i to the entropy density is then 
i 



g*s 



w 9l 



T / 4tt 4 



(3 
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e u ±l 



du 



e u ± 1 



-du 



(5) 



These integrals can be evaluated numerically once one 
knows which particle species should be included in the 
sums, and what value to take for the transition temper- 
ature T c . 

To answer the first question, consider the thermaliza- 
tion of the p(750) meson (in units of % = c = k = 1). 
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The /o(750) number density at temperature T = 0.2 GeV 
is 

N p = n spm {2TTm p T) 3/2 e~ m " /T w 0.065 GeV 3 . (6) 

The cms velocity in elastic ir + p — > 7r + p collisions is 
Wcms ~ 0.15. Taking the cross-section to be about 30 mb 
or a « 7.7 x 10 7 GeV -2 we find the reaction rate 

T = (N p v cms <j) » 7.5 x 10 5 GeV . (7) 

This should be compared to the decay rate of the rho 
meson, T p = 0.15 GeV. Thus there is plenty of time for 
rho mesons in the tail of their thermal distribution to 
thermalize in collisions against pions. 

One comes to the same conclusion when one studies 
the reaction 7 + p — > 7 + p : the reaction rate is then 
1.7 x 10 4 GeV for an approximate cross-section of 0.1 mb. 

One also comes to the same conclusion when one turns 
to heavier meson resonances and to baryon resonances. 
We therefore include all known meson and baryon reso- 
nances @ up to 3 GeV mass in the sums in g* and g*s 
at temperatures below T c . This is the procedure followed 
also by the earlier studies [2, La 0] ■ 

In spite of the strong exponential decrease in the num- 
ber density of non-relativistic species with heavy masses 
rrii > T c , such species will contribute significantly be- 
low T c . For example the proton and neutron contribute 
0.25 each to g*(T) at T — 150 MeV (versus an unatten- 
uated 3.5) even though their mass of roughly 939 MeV 
is 6.3 times higher. The sum of unattenuated mesonic 
and baryonic degrees of freedom below 3 GeV mass is 
gmesons = 495 and gbaryons = 320, respectively. As a 
result g*(T) rises quite rapidly for temperatures above 
T = 100 MeV. 

Near the transition temperature T c the ideal gas ap- 
proximation breaks down and the transition from the de- 
confined quark phase to the confined hadronic phase can 
only be studied with finite temperature lattice simula- 
tions. No reliable value for T c has yet been obtained, but 
it is believed to be near 170 MeV ||. Srednicki et al. 
proposed that the transition my be smooth and bracketed 
by the cases T c = 150 MeV and 400 MeV. Their proce- 
dure to obtain smooth curves in the context of the 'bag' 
model has only been elucidated by Sarkar in a footnote 
on p. 1513 |4| which refers to a private communication 
from K. A. Olive. But whether the change is smooth or 
discontinuous still remains open 

In Fig. ^ we plot g*{T) and g*s(T) as (almost) discon- 
tinuous functions of logT for the range 100 GeV > T > 
10 keV. For purposes of the figure we have somewhat ar- 
bitrarily chosen the discontinuity at T c = 200 MeV (as in 



3|) and a smoothing naive transition function that mixes 
the quark and hadronic phases over a range of roughly 
5 MeV on either side of T c . Since it is quite laborious 
to calculate the curves in Fig. ^ and difficult to read 
them precisely there, we give here some numerical ap- 
proximations over the range T = 155 — 300 MeV. For the 
hadronic phase we have 

g*(T) = 2.388 x 10" 4 T 2 ' 287 , max error 5.3% , 
g* s (T) = 3.631 x 10~ 4 T 2192 , max error 5.5% . (8) 
For the quark phase we have 

5»(T) = 24.604 x T 01713 , max error 0.2% , 

S*s(T) = 22.961 x T 01810 , max error 0.2% . (9) 

It is instructive to exhibit g*(T) in the quark phase 
and the hadronic phase separately, as we have done in 
Fig. H This shows a very well defined cross-over at 235 
MeV. Taking this temperature to be T c one obtains a 
curve which is indistinguishable from the curves in Fig.^ 
(because of the logarithmic scales). The cross-over tem- 
perature is defined quite precisely due to the steepness 
of the hadronic phase, but it does of course depend to 
some extent on the precision of the hadronic resonance 
energies, as well as on unknown particles and interactions 
above 1 TeV. 

Our most interesting conclusion is that a plain count 
of effective degrees of freedom sets an upper limit to the 
temperature of the quark-hadron transition. From the 
g*(T) cross-over one has 

T c < 235 MeV , (10) 
and from the g*s(T) cross-over 

T c < 245 MeV . (11) 

Higher T c values would give an unphysical spike and 
lower values would increase the discontinuous drop in 
g*(T) just below T c . 
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FIG. 1: Effective degrees of freedom for the energy density 
g*(T) and for the entropy density g*s{T). The quark-hadron 
transition is chosen to occur at 200 MeV. 




FIG. 2: Effective degrees of freedom g*{T) due to the quark 
and hadronic phases, respectively. For the quark phase we 
have taken T c = 150 MeV, for the hadronic phase T c — 300 
MeV. The cross-over occurs at 235 MeV 



